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Abstract. The article surveys quantization schemes for metric graphs with 
spin. Typically quantum graphs are denned with the Laplace or Schrodingcr 
operator which describe particles whose intrinsic angular momentum (spin) 
is zero. However, in many applications, for example modeling an electron 
(which has spin-1/2) on a network of thin wires, it is necessary to consider 
operators which allow spin-orbit interaction. The article presents a review 
of quantization schemes for graphs with three such Hamiltonian operators, 
the Dirac, Pauli and Rashba Hamiltonians. Comparing results for the trace 
formula, spectral statistics and spin-orbit localization on quantum graphs with 
spin Hamiltonians. 



1. Introduction 

Quantum graphs have become a widely applied model in mathematical physics. 
They fill a void between the simple quantum mechanics of the interval and more 
complex problems where it is often desirable to develop intuition in this quasi- 
one-dimensional system. Motivation for this approach comes from the picture of 
quantum mechanics on a three dimensional network of thin wires. It is expected 
that properties of the quantum graph will also be manifest in truly three dimen- 
sional systems, a subject of current research see e.g. |16l I19L 129] . Most work on 
quantum graphs has concentrated on the properties of the Laplace and Schrodinger 
operators, see |28j for a review. However, experiments on mesoscopic systems, sys- 
tems with lengths on the scale of nanometers, measure features of electrons confined 
to these narrow wires. Indeed quantum graphs were first proposed as models of 
free p-electrons in organic molecules like naphthalene, see |27j . The electron is a 
fermion possessing an additional quantum mechanical degree of freedom, intrinsic 
angular momentum called spin. The effect of spin in the particles physics is highly 
nontrivial - the connection between spin and particle statistics requires fermion 
eigenf unctions to be antisymmetric under particle exchange generating a wealth of 
phenomena in physics and chemistry. For single particle quantum mechanics the 
effect of coupling between the spin dynamics and particle dynamics - spin orbit 
coupling - also generates new physics observable in experiments. Understanding 
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how the presence of spin influences the quantum mechanics of graphs is therefore 
an important current question. 

In this article I review the present state of knowledge of quantum mechanics 
on graphs where the Hamiltonian operator acts on multi-component wave func- 
tions that describe particles with spin. The results fall into two broad categories, 
stemming from problems in quantum chaos and mesoscopic models of localization. 
A good introduction to the application of graph models in quantum chaology is 
found in |18j . One characteristic quantum phenomena associated with chaotic 
classical dynamics is an energy level spectrum with statistics that resemble those 
of an ensemble of random matrices, this correspondence is known as the Bohigas- 
Gianoni-Schmidt conjecture [7]. The choice of random matrix ensemble depends 
on the symmetries of the quantum system. For systems with half-integer spin time- 
reversal symmetry, in particular, takes a different form to that found in systems 
with integer spin. This results in a doubly degenerate spectrum, Kramers ' degen- 
eracy, and increased level repulsion between distinct eigenvalues. A notable success 
in the application of quantum graphs with spin to quantum chaos has been to ex- 
plain the effect of this change in symmetry in the spectral statistics of graphs with 
randomly chosen spin transformations at the vertices in terms of periodic orbits 
[111 110] . These arguments have recently been extended to transport properties of 
chaotic cavities with spin-orbit interaction |14j . 

The Rashba Hamiltonian was introduced to model the two-dimensional electron 
gas confined on the surface of a narrow gap semiconductor [34] . The effect of the 
Rashba spin-orbit interaction on a graph can be detected in localization phenomena. 
Numerical results for one-dimensional diamond chains show ideal localization, zero 
conductance, in the absence of a magnetic field [2j [3]. For the two-dimensional 7jj 
lattice localization is not perfect but the spin-orbit coupling and magnetic field can 
be tuned to suppress the conductance. Analytical results for the spectrum show 
Bloch bands that collapse to infinitely degenerate eigenvalues [33) . 

The article is organized as follows: Section [2] introduces the nomenclature of 
quantum graphs. In section [3] we collect results on the quantization of graphs with 
spin Hamiltonians. Section [4] discusses the trace formula for a quantum graph with 
spin and in the balance we describe effects of spin-orbit interaction that have been 
investigated using quantum graphs in spectral statistics, section [5l and Rashba 
localization, section [6l 



2. Quantum graphs 

I have endeavored to follow standard terminology for quantum and metric 
graphs throughout. A reader already familiar with these ideas is encouraged to 
skipped ahead and refer back to this section if necessary. 

A metric graph G consists of a set V of vertices and a set £ of edges. N = \£\ 
denotes the number of edges. Each edge e = {v, w} connects a pair of vertices 
v, w £ V and the edge e is associated with an interval [0, L e ] where L e e [L m - lu , L max ] 
is the length of the edge, < L m - m < L max < oo. Introducing a coordinate 
x e G [0, L e ] on the edge the graph naturally becomes directed with the initial 
vertex v of the edge located at x e = and the terminal vertex w at x e = L e . A 
directed edge will be denoted with parentheses, e = (v, w). Functions on the graph 
are defined by the set of functions on the edges / = (/*, . . . , / ). 
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The quantization of the graph requires the introduction of the Hilbert space, 
(2.1) L 2 (G)=0L 2 (O,L e ). 

Other spaces of functions are defined analogously, for example W 2,2 (G) is the space 
of functions on G with components drawn from the L 2 -Sobolev spaces W 2,2 (0, L e ). 
In this article we will be concerned with the properties of operators acting on multi- 
component functions on G where the relevant Hilbert space is then L 2 (G) (g> C n for 
spinor valued functions with n components. The scalar product of /, g is defined 
by the scalar product on the edges, 

(2-2) (/,5>=E</ e >2 e >e 

e££ 

where (f e ,g e ) e is the scalar product in L 2 (0,L e ) ® C n . 

3. Quantization schemes 

The approach taken to quantizing a graph follows the same pattern for all the 
spin operators. One takes the relevant operator in three dimensions, restricting it to 
one dimension to define the operator on an interval. As spin is an intrinsic angular 
momentum naturally connected to rotations in three dimensions it is possible to 
restrict such three dimensional operators in different ways. However, if the lengths 
of the edges can be chosen freely it is always possible to embed a graph in three 
dimensions in which case the three dimensional notion of a particles spin still makes 
sense. This is also the natural physical situation to study where a quantum graphs 
is used to model networks of thin wires. Having defined the operator on the edges 
matching conditions for the multi-component wave functions at the graph vertices 
are specified to make the operator self-adjoint. 

3.1. The Dirac operator. The free Dirac operator on a graph was first con- 
sidered by Bulla and Trenkler |15] . They take a two dimensional Dirac operator on 
the intervals, construct self-adjoint realizations and analyze certain special cases. 
In [9] we introduced an approach which extends the classification of self-adjoint 
Laplace operators |24) to classify self-adjoint realizations of the Dirac operator in 
terms of matching conditions for both two a four component spinor wavefunctions. 
Further quantization schemes for a metric graphs with a Dirac operator have been 
developed by Post. In 32 he introduces the notion of a first order supersymmetric 
Dirac operator on both discrete and metric graphs in a very general framework and 
calculates the index of the operator proving that it agrees with the corresponding 
index of the discrete Dirac operator. First order boundary triples are used to define 
a Dirac operator in [311 . In the following, for simplicity, we follow the approach we 
developed in [9]. 

In one spatial dimension the Dirac equation is 

d ( d \ 

(3.1) ih—&(x,t)= -iHca — + mc 2 (3 )&(x,t) , 

at \ ox J 

where a and (3 are matrices that satisfy the Dirac algebra a 2 = 1 = I and a (3 + 
(3a = 0. It is natural to consider this equation on the line as deriving from a 
restriction of the Dirac operator in three spacial dimensions with the consequent 
notions of particles, antiparticles and spin all being carried over. The equation 
then arises from implementing the Poincare space-time symmetries in relativistic 
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quantum mechanics, see 
convenient choice being 



and the matrices a and /3 are four dimensional, a 



(3.2) 



(3 = 
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Alternatively one might look for the simplest faithful irreducible representation of 
the Dirac algebra which uses two dimensional matrices. Here we consider only 
the first approach with four component spinor wave functions. The quantization 
of two component wave functions was carried out following the same scheme in 
[9] , We found that to make such a Dirac operator time-reversal symmetric (in 
a nontrivial way) it is necessary to replace each edge of the combinatorial graph 
with two directed bonds one running in each direction, effectively introducing four- 
components to the wave function on each edge. The scattering matrices at the 
vertices and the spectrum of the two component operators on pairs of directed 
bonds are then equivalent to those of a four component Dirac operator. 
A free Dirac operator on an edge has the form 



(3.3) 



= —\hca 1- mc fi . 



T>° is defined as T> with the domain W ^(G) £g>C 4 , the Sobolev space of functions / 
on the graph where the components (and their first derivatives) are L 2 (0,L e ) and 
the functions vanish at graph vertices, / e (0) = f e (L e ) = 0. Self-adjoint realizations 
of T> are defined by extensions of the closed symmetric operator T>° to domains that 
are subspaces of Wq' 1 (G) £g> C 4 isotropic with respect to a skew Hermitian form 
n(f,g) = (Vf,g)-(f,Vg). 

Integrating by parts one may rewrite f2 as a complex symplectic form on C 8N 
depending only on values of the spinors at the graph vertices. 



(3.4) 



fi(/,s) = ( f 



h B 

-Lib 



where 

f+ = (A 1 (G), #(0), ■ ■ • , /f», jf(Q), fl(Li), /!(£i), . . . , f?(L N ), f 2 N (L N )) , 

f T = (-/i(0), ft(0), -/f (0), / 3 », /Kir), -fl{L,), f?(L N ),-f 3 N (L N )). 

The vector f+ contains the first and second components of the spinors and f_ 
the third and fourth components. Maximal isotropic subspaces with respect to f2 
correspond to maximal subspaces of vectors in C SN on which the symplectic form 
(|3.4[) vanishes. 

As in the case of the Laplace operator a linear subspace of vectors in C 8N will 
be defined as those vectors satisfying the equation 



(3.5) 



Af, 



Bf_ = , 



with complex matrices A and B. In order for the space to have dimension AN 
the 4iV x 8./V matrix (A, B) must have maximal rank. Kostrykin and Schrader 
|24) showed that such a subspace is maximally isotropic if and only if (A, B) has 
maximal rank and AB^ is Hermitian. 
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Theorem 3.1. Af+ + Bf_ = defines matching conditions of spinor wave 
functions at the vertices of a graph for a set of N bonds. These matching conditions 
yield a self-adjoint realization of the Dirac operator on the graph if and only if 

rank(A, B) = AN and AB f = BA f . 

When considering graphs it is often convenient to work with scattering matrices 
at the vertices rather than with the matching conditions directly. To obtain the 
vertex scattering matrix of the Dirac operator we note that eigenfunctions on an 
edge, the solutions of T>ip e (x e ) = Eip e (x e ), are plane waves. 

( 1 \ ( \ 





(3.6) 
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for k > and positive energy E and j(k) given by 

E — mc 2 



(3.7) 



7 (fc) = 



hck 



E = ^/(hck) 2 + m 2 c 4 



Consider a vertex v of valency d v whose matching conditions are defined by a pair 
of matrices AW,BW, so A^xp + + M^ip_ = 0. To simplify the formulae we take 
edges meeting at v to be aligned with v at x e = 0. We may define two vectors of 
plane-wave coefficients for the incoming and outgoing waves respectively, 



(3.8) 



7? 



(Aqj A/3; 



n d n d ) T 



u d , 



fi + fi and tp_ 



4) ■ 

vy(k)(~JjL — n). Consequently, 
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Then ip + 

(3.9) V 
The vertex transition matrix 

(3.10) = -(A^ - i 7 (fc)B^)- 1 (A^ ,) + i 7 



v) 



relates coefficients of incoming and outgoing plane- waves at v. It is unitary as 
j^MjjMt = B^^A^^, which ensures current conservation at the vertex. 

3.2. Time-reversal invariance for the Dirac operator. Time-reversal 
symmetry in a system with half-integer spin takes a different form from the sym- 
metry under complex conjugation familiar for scalar functions, see [20] . The time- 
reversal operator T for systems with half-integer spin remains anti-unitary but it 
squares to minus the identity, T 2 = —I. For any Hamiltonian H invariant under 
time reversal T commutes with H and if ip is an eigenfunction so is Tip. However, 
when T 2 = —I the two eigenfunctions ip and Tip are orthogonal and eigenvalues of 
H are doubly degenerate - known as Kramers' degeneracy. 

For the Dirac operator on a graph to be invariant under time-reversal the 
operator T must commute with the Hamiltonian which forces the mass term to 
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vanish fixing j(k) 
symmetry implies 

(3.11) 

where 

(3.12) 



1. For the vertex transition matrices we find time- reversal 



(T^f = -(I dv ® J)?W fti.® J) 



1 

-1 



If we consider the 2x2 block of TW that relates the pair of incoming spinors 
arriving from edge e to the outgoing pair on / p. lip is equivalent to, 

(3.13) (TW)/ e = \(T^ v) ) ef \ ((T^) 6 ^' 1 



This suggests a simple method of constructing time-reversal symmetric transition 
matrices that relates transition matrices of the Dirac operator to those of the 
Schrodinger operator on the graph. Let 

(3.14) T (u) = (U { " ) )~ 1 (XW ®Ia)U w , 

where U^) is a block diagonal matrix U^' = diagjui, . . . ,Ud v } with Uj € SU(2) 
and X is a unitary symmetric d v x d v matrix. Then respects p. lip and results 
from a self-adjoint realization of the Dirac operator on G. X can be thought of 
as the transition matrix of the graph with out spin. For the Laplace operator on 
the graph time-reversal invariant matching conditions at the vertices generate a 
symmetric unitary vertex transition matrix. The elements Uj G SU(2) generate 
spin rotations at v. We may now define, for example, Neumann like matching 
conditions for the Dirac operator at the vertex v. 

(uW)- 1 C? = («g))- 1 C?=£ f («) eue^v 

(3.15) 



4" ) )~ 1 c = o 



The sum is over edges e connected to v, e 
in the form (|3.5[) by matrices 
(3.16) 



v. These conditions can be represented 
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Substituting the matching conditions in (|3 . 1 0|) produces a transition matrix of the 
form (|3.14p with [X^^] e / = 2/vd — § e f the familiar transmission amplitudes of the 
Laplace operator on a graph with Neumann matching conditions at the vertices, 
see 



3.3. The Pauli operator. For comparison we introduce a quantization scheme 
for a simple form of the Pauli operator which is closely related to the Laplace quan- 
tization [24] . Consider the following operator on the edges of the graph, 

< 3 - 17 > v < - -s 



B, 
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where er is the vector of Pauli matrices and B e is a constant vector depending on 
the edge. The operator acts on two component functions f e in the Hilbert space 
L 2 (0,L e )®C 2 . 

V° is defined as V with the domain Wq ' 2 (G)(g>C 2 , the Sobolev space of functions 
/ on the graph where the components (and their first and second derivatives) are 
in L 2 (0,L e ) and / e (0) = f e (L e ) = f' e (0) = f e (L e ) = 0. Self-adjoint realizations of 
V are again defined by extensions of the closed symmetric operator V° to domains 
that are subspaces of W 2,2 (G) ® C 2 isotropic with respect to a skew Hcrmitian 
form n(f,g) = (Vf,g) - (f,Vg). As (B e .er)t = B e .cr the form Cl(f,g) is the same 
as that obtained for the Laplace operator where the number of edges is doubled. 
A domain on which the operator is self-adjoint is specified by a pair of 4N x AN 
matrices A and B, see |24) . where (A,B) has maximal rank and AB^ is Hermitian 
via the linear equation, 

(3.18) Af + Bf' = 0. 

f and f ' are vectors of the components of / and its derivative evaluated at the ends 
of the intervals. 

f = (/i(0), /» ■ • • , /!», / 2 », fl(L e ), fl{L e ), f» (L e ), f 2 N (L e )) T 
f = (tf(0), /^(0) • . • , f[ N (0), ff (0), -f x l {L e \ f?(L e ), . . . , -f[ N (L e ), ff(L e )f 

To define the vertex transmission matrix it is convenient to diagonalize the spin 
transformation on the edge, 

(3.19) B e .CT = IBelWeCTa^ 1 , 

where u e is in SU(2). Then plane wave solutions for the eigenproblem V e ip e = 
are given by 

(3 20) f -u ( ^1 + Afe ~ ifc e 

(fc|) 2 = {kf) 2 + 2|B e | and A = (k\) 2 + |B e |. To obtain the vertex scattering 
matrix consider a vertex v of valency d v with the edges aligned so v is the initial 
vertex of each edge. Vectors of incoming and out going plane-wave coefficients are 
respectively, 

*i* = (Ai»A2>--->Ai,A2) T , 

Then defining K (v) = diag{fcj, fcg, /c 2 , . . . , /c^"} we nave 

(3 22) ^ = C/W(7? + V), 

ip' = -iuMK<- v \~jt -*jt) , 

where LA") is again a block diagonal matrix, 17 W = diag{ui, . . . , Ud v }- Conse- 
quently substituting in (|3.18p the vertex scattering matrix is 

(3.23) = -(U^y 1 (A^ - iK^M^Y (a (v) + iK^M^A U {v) . 

In this form the transition matrix has built in spin rotations at the vertices of a 
similar form to that described in (|3.14p for the Dirac operator. A Pauli operator 
can also be defined for systems with higher spins, acting on multi-component wave 
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functions of higher dimensions. This was used in [llj to extend results on the 
spectral statistics of the Dirac operator to any spin. 

3.4. The Rashba Hamiltonian. The Rashba Hamiltonian, introduced in 
[34j to model electronic structure of confined narrow-gap semiconductors, has also 
been studied as an operator on quantum graphs. Harmer [22 constructs the scat- 
tering matrix of the Rashba Hamiltonian on a ring with an arbitrary number of 
semi-infinite wires attached to investigate spin filtering. In section [6] we will dis- 
cuss localization phenomena generated by Rashba spin-orbit coupling [2j [3], 133) . 
So to complete our survey of quantization schemes we now turn to the Rashba 
Hamiltonian on a graph. 

With the Rashba Hamiltonian it is useful to have in mind a graph with straight 
edges embedded in the x-y plane in R 3 . We will denote by et vw -s — (ef vw \ x , e/ vw \ y , 0) 
a unit vector in the direction (vw). Let Ar vw \ denote the magnetic potential on 
(vw), the component of a magnetic vector potential in the direction e^). For 
simplicity we take A( vw ^ to be constant, as is the case with a uniform magnetic 
field B where Ar vw j = ^(B x v).e( vw ) and v is a vector locating the vertex v. 
is the Rashba coupling constant and n = (0,0, 1) is a unit vector perpendicular to 
the x-y plane. The Rashba Hamiltonian on an edge e reads, 



A self-adjoint realization of the operator with delta type matching conditions at 
the vertices was defined by Pankrashkin in [33] . The matching conditions at the 
vertices are analogous to delta type conditions of the Schrodinger operator where 
e(v) — defines Neumann like coupling at the vertex v. 

Theorem 3.2. Denote by an operator 1Z the Rashba Hamiltonian on the edges 
of the graph acting on functions f £ W ' (G) (g) C 2 satisfying at each vertex v G V, 



(3.24) 




- A - ( CT x n ) e » 



cr denotes the vector of 2 x 2 Pauli matrices and 



(3.25) 




f(«w)(0) = f(uv){L( uv )) = f(w), (vw), (uv) e £ , 




i(A uv) i ^'R^(uv) 

) WO) 

(uv) J 




Then 1Z is self-adjoint. 



4. The trace formula 



A trace formula for the Laplacian on a graph was first produced by Roth 
[35] . Here the trace formula for the density of states of the free Dirac and Pauli 
operators is stated in the form developed by Kottos and Smilansky |25L 126] for the 
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Schrodinger operator. A full review of the state of the art for the trace formula of 
a Schrodinger operator on a graph can be found in this volume [8 . 

The spectrum of a quantum graph can be computed using the scattering matrix 
of the graph. The scattering matrix is a unitary AN x AN quantum evolution 
operator acting on the vector of all plane- wave coefficients. It is defined via 

(4.1) (fc) = S vw a {uv)(vx) u (uv)(vx) e ifci <-> , 

where a {uv){vx) = |( T M)(™)(-)| and u {uv){vx) = (iW) (T W ) (m)(ra) | . 
The fc-spectrum of the Dirac operator are the solutions of the secular equation, 

(4.2) |I-S(Jb)|=0. 

This quantization condition is obtained by considering vectors of plane-wave coef- 
ficients invariant under multiplication by S(fc) which define eigenfunctions on the 
graph. 

The density of states is a comb of delta functions located at wave numbers fc„, 
d(k) — YlnLi ~ The trace formula relates this spectral distribution to the 
set of periodic orbits on the graph. A periodic orbit p of n steps is a sequence 
p = (ei, 62, . . . , e„) of directed edges on the graph, where the terminal vertex of ej 
is the initial vertex of ej+i and e„ terminates at the origin of e%. We will denote 
the set of Periodic orbits of n bonds P n and and l p is the metric length of the 
orbit p. The total length of the graph is L = J2 e es ^e- A periodic orbit p may be 
a repartition of a shorter primitive orbit in which case r p denotes the number of 
repartitions of the primitive orbit used to produce p. 

Theorem 4.1. For the Dirac operator on a graph the density of states is, 

\ i 

d(k) = 1 — > — A p tr(d p ) cos(fcL) , 



where 



7T " r n 

p 



®~ e\e2®~ eves • ■ 



' U 



The trace formula has the same form as that of the Schrodinger operator on 
a graph. The first term is a Weyl term, the mean density of states. The second 
contribution, the oscillating part, contains information about correlations in the 
spectrum written as a sum over all periodic orbits. The difference from the scalar 
version appears in the additional trace of d p e SU(2) which is the product of 
the spin transformations at the vertices picked up while traversing the orbit p. A 
similar term appears in the semiclassical Gutzwiller type trace formula for the Dirac 
operator [HI [13]. The trace formula is derived from the secular equation (|4.2p in 
[9|. For the Pauli operator on the graph a similar formula holds and the elements 
ybjbj+i g SU(2) may be replaced by a higher dimensional irreducible representation, 
see [TT1. 



5. Spectral statistics 

Quantum graphs with spin Hamiltonians provide a simple model in which to 
investigate the effects of spin-orbit interaction. Our first example is the effect of 
spin-orbit interaction on the statistics of the graph spectrum. 
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Figure 1. The integrated nearest neighbor spacing distribution 
and a histogram of the spacing distribution for a Dirac operator 
with Neumann matching conditions on a fully connected graph 
with four vertices, 24234 levels. Figure taken from [9]. 



For systems with half-integer spin and time-reversal symmetry eigenvalues of 
the Hamiltonian are doubly degenerate. When discussing spectral statistics it is 
therefore convenient to lift Kramers' degeneracy and also rescale the spectrum so 
that the mean spacing is one, for more details see |20j . In the following we will 
assume such an unfolding of the spectrum has already taken place. 

The spectrum of the Dirac operator can be investigated numerically using the 
secular equation (|4.2[) . A readily available statistic for computation is the nearest 
neighbor spacing distribution p(s), the probability density of the spacing s between 
consecutive levels. Figure [T] shows a histogram of the nearest neighbor spacing 
distribution and the integrated distribution for the Dirac operator on a fully con- 
nected square. The spin transformations at the vertices were chosen randomly from 
SU(2). The numerical results show good agreement with the spectral statistics of 
the Gaussian symplectic ensemble (GSE) of random matrices - the ensemble of 
Hermitian matrices with probability measure invariant under symplectic transfor- 
mations. This correspondence with random matrix statistics, seen in a wide variety 
of chaotic quantum systems, is known as the Bohigas-Giannoni-Schmidt conjecture 
7J . The point of particular interest to us here is the ensemble of random matrices to 
which the statistics of the quantum spectrum correspond. For systems with time- 
reversal symmetry and half-integer spin like the Dirac operator on a graph one 
generically sees GSE statistics, while for integer spin the corresponding random 
matrix statistics are those of the Gaussian orthogonal ensemble (GOE). 

The two point correlation function of the unfolded spectrum is, 

(5.1) R 2 {x)=\ira\j d{y) d(y - x)dy - 1 . 

where d(y) is the unfolded density of states. The fourier transform of i?2, the form 
factor, has been the spectral statistic most amenable to theoretical investigation. 
Using the trace formula for the Dirac operator (|4. 1[) the spectral form factor can 
be written as a sum over pairs of periodic orbits, 

(5.2) KmrM = i ^r 5 Yl 77 e^ + ^ tr(dp) tv(d q ) 8 (r - ±\ 5 lp , lq 
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for r positive [9] . This is closely related to the form factor of the Laplace operator 



(5.3) K^Ur) - £ ^A P A, S (r - i) . 

Generically, for time-reversal symmetric boundary conditions, the form factor of 
the Dirac operator converges to the limiting distribution of the GSE when one 
considers the limit of a sequence of graphs with increasing numbers of edges. In 
the same way the form factor of the Laplace operator on the graph converges to 
the GOE form factor. Comparing power series expansions of the random matrix 
form factors for t < 1 shows that they are also closely related 



^gse(t) 

(5-4) 



Calling K m the term containing r m the relationship can be written 

j \ m+1 



T 


T 2 


r 3 


r 4 


2 + 

T 


T 2 


r 3 


+ 12 

r 4 


2 ~ 


4 




12 



(5-5) K% se (t)={--) K% oe (t). 

In (|5.2[) and (|5.3[) a pair of orbits p, q only contributes to the form factor if 
lp = l q . Much of the success studying the form factor via periodic orbits has 
come by identifying sets of orbit pairs where the partner orbit of p is generated by 
permuting the arcs of p between self intersections and or reversing the direction of 
arcs. Clearly permuting the order of edges of an orbit on a graph or reversing the 
direction a sequence of edges is traversed does not change the orbit length and pairs 
constructed in this manner must contribute to the form factor. The first example of 
this approach in a chaotic quantum system was the diagonal approximation used by 
Berry [6] to obtain the first term in the power series expansion of the form factor by 
evaluating the contributions of an orbit with itself, p — q or with its time reversed 
partner using the sum rule of Hannay and Ozorio de Almeida [21] . Sieber [36] and 
Sieber and Richter |37| were the first to evaluate the second order contributions by 
considering orbits with one self-intersections. This was extended to higher orders 
for quantum graphs in [4l [5] and recently for general contributions of orbit pairs 
with any number of self-intersections [30 . 

To evaluate such contributions to the form factor of a quantum graph requires 
us to consider the combined limit L — ► oo, r — > with Lt — > oo while the mean bond 
length L = L/N is kept constant. In this limit long orbits dominate the sum and 
from these the proportion with r p ^ 1 tends to zero and may safely be excluded. 
We also approximate l p w nL where n is the number of bonds in p. To isolate 
a spin contribution to the form factor (|5.2p we assign the spin transformations 
iO randomly with Haar measure from a subgroup r C SU(2) so that the only 
correlation between the stability factor A p and the spin rotation d p comes from the 
common structure of the orbit. Let D be a set of pairs of orbits (p, q) of length 
n where each partner orbit q is obtained from p by the same permutation and or 
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reversal of arcs of p. We can factor the sum over the orbit pairs in D so, 
rL\2 Ap A q ^ + ^hv(d p )tv(d g ) 



(5.6) 



4 E A p A i ^ + ^) (r^r J2 tr W tr(d q ) 



\D ^ " « / \ D 

1 1 D / V ' D 



For long orbits and randomly chosen spin transformations we argue that 
(5.7) — E tr(dp) tr(d,) = (tr(dp) tr(d g )) r 

' ' D 

the average over r with Haar measure. This provides a sketch of how a spin 
contribution for certain sets of orbits can be isolated from terms in the form factor 
sum. The remaining term in (|5.6|) is the contribution of these orbit pairs to the 
form factor of the Laplace operator on the graph. 

In [101 111] we calculate the spin contribution of such sets of orbit pairs. 

Theorem 5.1. Consider a pair of orbits p,q on a graph where q is related to 
p by a permutation of arcs of p and or reversing the direction of arcs. Let the spin 
transformations on the graph generate a subgroup r C SU(2). 



(dp : dqj 



where the average is with respect to Haar measure. The exponent m — 1 is the 
number of self intersections at which the orbit p has been rearranged to produce q. 

Multiplying by the factor of 1 /4 from the form factor of the Dirac operator we 
obtain (— ^) m+ ■ This generates the same relationship as exists between the form 
factors of the GSE and GOE (|5.5p provided that pairs of orbits where the partner 
q differs from p at m — 1 self intersections contribute at order r m in a power series 
expansion of K(t) . This is precisely the relationship found in derivations of the form 
factor without spin from periodic orbit theory [H [5J [fil I30|, 1361 137] . The spectral 
statistics of quantum graphs with spin are investigated using a completely different 
technique in [17) . Here the authors compute spectral correlations using super 
symmetric arguments and obtain the characteristic GSE statistics for a quantum 
graph with spin. However, as the model of spin dynamics is specific it can't suggest 
conditions on the spin model necessary for Random matrix statistics. 

6. Rashba localization 

The effect of a spin orbit interaction can be found not only in spectral statis- 
tics but also in macroscopic properties of quantum systems. Spin orbit coupling 
in quantum graphs can be used to generate localization phenomena in periodic 
structures. In [2j [3] Bercioux et. al. demonstrate localization, zero conductance, 
numerically for a diamond chain lattice with equal edge lengths, figure [5J They 
consider the Rashba Hamiltonian on the graph where the chain is embedded in 
the x-y plane and the magnetic field is perpendicular to the plane, as described in 
section l3Tl In the absence of a magnetic field the Rashba coupling can be tuned to 
so that spin transformations along the graph edges conspire to produce cancelation 
in the wave functions at the vertices of valency four. An electron traveling across 
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Figure 2. A schematic view of the diamond chain. 

(a) » 

1.5 

(G) ki : 

0.5 
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2k R Lir- 1 2k R LTT~ 1 

Figure 3. (a) Conductance averaged over k- m as a function of 
the spin-orbit coupling strength for the diamond chain (continu- 
ous line) and for the ladder (dashed line). The two finite-size sys- 
tems connected to input/output leads are shown in the inset. The 
parameters used for the calculation are: 50 elementary loops, k- ln 
uniformly distributed in [0, 7r/L]. (b) Conductance as a function of 
the spin-orbit coupling strength for the diamond chain (continuous 
line) and for the ladder (dashed line) for a fixed value of fc; n = fcp . 
The parameters used for the calculation are: 50 elementary loops, 
kpL = nir + 2. being n an integer. Figure courtesy of Bercioux et. 
al. ®. 




one of the cells, from v to w in figure[2j can take two possible paths. Along an edge 
e it acquires a spin rotation 

(6.1) d e = exp [ — \2k R L e {cr x n).e] . 

along with the dynamical phase e lkL ". One may consider the spin transformation 
(|6.ip as a non-Abelian phase which gives rise to destructive interference only for 
the special case of the square. Although the effect only appears for this particular 
geometry the phenomena is seen numerically to be stable under small perturbations 
of the edge lengths about a mean where they are all equal. Consequently it would 
be amenable to experimental investigation [2]. 

Figure [3] (a) shows the conductance across a diamond chain and a ladder as a 
function of the spin-orbit coupling averaged over the injection energies. Localization 
in the form of destructive interference at the valency four vertices can also be 
generated by the Aharonov-Bohm phase in the absence of spin-orbit coupling k R = 
for particular values of the magnetic field strength. These vertices which bound 
an initially localized wave-packet are known as an Aharonov-Bohm cage. The same 
cage generated via the Rashba effect by tuning the spin-orbit coupling in the absence 
of a magnetic field. 
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(a) (b) 2kuL-K 1 (spin-orbit coupling) 




4>/4>Q reduced flux 2knLn 



Figure 4. (a): A finite-size piece of the 7jj network connected to 
reservoirs. There are 4 input and 4 output single-channel leads. 
The number of input channels, considering spin, is N{ n = 8. (b): 
Averaged conductance per channel, (G)k/Ni n , as a function of the 
reduced flux evaluated with kji = (solid line), and of spin-orbit 
coupling evaluated with 4>/c/)q — (dashed line) for the T3 lattice 
with 200 quantum wires (89 rhombi). (c): Averaged conductance 
per channel, (G)fc/Mn, as a function of the reduced flux evaluated 
at 2knLpi~ 1 = 0.5 (solid line), and of spin-orbit coupling evaluated 
at (f>/(j)o = 0.5 (dashed line), (d): Grey-scale plot of the averaged 
conductance per channel, (G)k/N- m , as a function of the reduced 
flux and spin-orbit coupling. Figure courtesy of Bercioux et. al. 
®. 



The spectrum and conductance of a Rashba Hamiltonian on the T3 lattice, see 
figuredfa), was investigated numerically in [3]. Spin-orbit coupling in this geometry 
no longer results in perfect localization but the Rashba effect can be tuned with 
or without the magnetic field to suppress the conductance, see figure |4jc) . Perfect 
localization is still possible when fcp = 0, figure 0Jb). 

In [33 the spectrum of the infinite T 3 lattice was characterized analytically via 
Bloch theory. To give a flavor of the result we restrict to the case of the free Rashba 
Hamiltonian. The matching conditions at vertices characterized in theorem 13. 21 are 
restricted to at most two types with a coupling constant e(v) — A at vertices of 
degree 6 and s(v) = fi at vertices of degree 3. An interesting situation occurs when 
cos kfi = 0, i.e kR S n/2 + 7rZ and the magnetic flux through a primitive rhombus 
uj = — 71/6 + 7rZ. The spectrum may be characterized as the set of values of E e R 
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such that 







cos \ E 



k 2 R 



6VE 



sin \J E + kj^j . ^cos \J 



0. 



6 



E + kl 



U 



3VE 
1 2 
2' 3 



sinA/S + ki 



along with solutions of the equation sin \JE + k R = 0. In this situation together 
with the expected bands there appear infinitely degenerate eigenvalues. Restricting 
further to Neumann like matching conditions, i.e. A = \i = the spectrum is the 
set of values of E for which, 



cos 2 ^E + k R G 



0,- 



u 



u 



1 2 

2' 3 



U{1}. 



Solutions of the equation 



cos \l E + k 2 R 



cos \l E + k R 



also appear as eigenvalues in the case of perfect localization induced by the magnetic 
field when u G it/2 + ttZ and kn G Z. Then all the bands collapse to a sequence of 
infinitely degenerate eigenvalues. It is notable that by tackling localization problems 
on graphs the simplicity and symmetry of the problems can be exploited to achieve 
elegant analytic results. 



Final remarks 

The article has been intended to survey techniques used to quantize graphs with 
spin. Our motivation goes back to the origin of quantum graph models. They were 
first proposed to describe electrons in organic molecules and modern experiments 
on mesoscopic networks also deal with the properties of electrons. Our examples of 
applications of spin Hamiltonians show that the presence of a spin-orbit interaction 
on graphs has a significant influence on the models features. A different class of 
spectral statistics or new modes of localization sensitive to the geometry of the 
network are both explained in terms of the spin rotations around periodic orbits. 

To take graph models with spin further there remain important open questions. 
Much has already been discovered about the convergence of Laplace operators on 
thin networks to the operator on a quantum graph and these results are sensitive to 
the boundary conditions at the surface of the domain. Quantum graphs with spin- 
orbit interaction should provided more interesting models of mesoscopic systems. 
To use them in this context it will be important to know when spin Hamiltonians 
converge to operators on quantum graphs. It will also be interesting to discover the 
physically relevant matching conditions for the various spin Hamiltonians. Then 
the ideal choice of model should introduce spin-orbit interactions in a way that 
maintains the simplicity that drew people to study graph models and provides real- 
istic results that can guide future analysis and mesoscopic experiments. So far the 
models generally fall into two classes: those where the spin transformation is at the 
vertices specified by matching conditions or those where spin rotation occurs along 
the edges. It seems likely that these may be equivalent as the convergence results 
for operators on thin networks will likely produce matching conditions dependent 
on the graph geometry at the vertices. 
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